Enhancement of Blackbody Friction due to the Finite Lifetime of Atomic Levels 



G. Lach, 1 M. DeKieviet, 1 and U. D. Jentschura 2 

1 Physikalisches Institut der Universitdt, Albert- Ueberle-Strasse 3-5, 69120 Heidelberg, Germany 
2 Department of Physics, Missouri University of Science and Technology, Rolla, Missouri 65409-0640, USA 

The thermal friction force acting on an atom moving relative to a thermal photon bath is known to 
be proportional to an integral over the imaginary part of the frequency-dependent atomic (dipole) 
polarizability. Using a numerical approach, we find that blackbody friction on atoms either in 
dilute environments or in hot ovens is larger than previously thought by orders of magnitude. This 
enhancement is due to far off-resonant driving of transitions by low-frequency thermal radiation. 
At typical temperatures, the blackbody radiation maximum lies far below the atomic transition 
wavelengths. Surprisingly, due to the finite lifetime of atomic levels, which gives rise to Lorentzian 
line profiles, far off-resonant excitation leads to the dominant contribution to the blackbody friction. 
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Introduction. — In Ref. PQ , the thermal drag force on an 
atom moving through a thermal bath at velocity v has 
been calculated on the basis of the fluctuation-dissipation 
theorem. In a nutshell, the fluctuation-dissipation the- 
orem states that any thermal fluctuation of a physical 
quantity (say, the electric field at finite temperature) is 
accompanied by corresponding fluctuations in the con- 
jugate variable (here, the atomic dipole moment) pro- 
vided the susceptibility (in the current case, the atomic 
polarizability) has a nonvanishing imaginary part. The 
imaginary part describes a dissipative process, in which 
the atom absorbs, then spontaneously emits, electromag- 
netic radiation. The dissipative fluctuations give rise to 
a drag force calculated using the Green-Kubo formula, 
as thoroughly explained in Ref. pQ. Further physical in- 
sight can be gained if one understands the process in 
terms of the direction-dependent Doppler effect [2]. The 
atom absorbs blue-shifted blackbody photons coming in 
from the front, while emitting these photons in all direc- 
tions, thereby losing kinetic energy due to net drain on 
its energy, and, as a consequence, on its momentum [3J. 

Here, we show that even more intriguing problems arise 
when one tries to evaluate the effect numerically, for sim- 
ple atoms. In atomic physics, the width T n for each indi- 
vidual energy level n needs to be determined separately. 
Atomic transitions can be driven even very far from res- 
onance, albeit with small transition probabilities, The 
blackbody spectrum is distributed over the entire fre- 
quency interval uj £ [0, oo), which leads to significant 
non-resonant contributions to the thermal friction. 

The authors of Ref. p] use correlation functions for the 
thermal electromagnetic fluctuations [H |S], in order to 
calculate the friction force acting on neutral, polarizable 
objects moving through uniform and isotropic thermal 
radiation. According to Eq. (12) of Ref. [1], the effective 
friction (EF) force, which acts in a direction opposite to 
the velocity v, is given as a spectral integral, 
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where /3 — l/(fcsT) is the Boltzmann factor and a{ui) is 
the dynamic polarizability of the atom. We here argue 
that the inclusion of the resonance widths due to the 
finite lifetimes of atomic levels is crucial in the calculation 
of the friction force. SI mksA units are used throughout 
this work. 

Narrow and finite width. — If we assume that all atomic 
transitions are infinitely narrow (of width e), then 
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where fo n denotes the oscillator strength of the transition 
and ojQn is the angular frequency for the transition from 
the ground state |0) to the excited state |n). In view of 
the Dirac prescription l/(x — ie) = P(l/x) + inS(x), the 
imaginary part of the polarizability is approximated as a 
sum of Dirac S peaks, 
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However, if one includes the width T n of the excited 
states, then the starting expression (see Chap. 8 of 
Ref. [B]) for the dynamic polarizability reads 
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Here, the decay width T n (uj) may be a function of the 
driving frequency uj. In a number of places in the lit- 
erature [e.g., see the text after Eq. (2) of Ref. |7J], it is 
assumed that 



r„(w) = T n (u) 
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One can justify the ansatz ^ in two ways (i) and(ii). 
(i) One may invoke an analogy with a damped, driven 
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harmonic oscillator, whose Green function g(t — t') fulfills 
the defining differential equation 
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so that the Fourier transform of the Green function reads 



as g{u) = 1 /(lo 2 — lo 2 — ijui), with 
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Assuming that T n (uj) = r„(w), this is proportional to 
the expression in Q under the obvious identification 
loq —} u>f) n , 7 — > r„. (ii) The decay width T n enters 
the propagator denominators in Eq. Q by a summation 
of self-energy insertions (5). The imaginary of the self- 
energy, divided by h, equals the decay width [91. The 
velocity-gauge expression QUI QT] for the decay rate, 
at resonance lo = io n and off resonance (for general lo), 
reads as 
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where p is the momentum operator, and and ^ n are 
the wave functions of the ground and excited state. 

By contrast, the so-called length gauge expression [8, 
fTOl ITT] for the decay width off resonance reads as 
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For atoms, using the commutator relation p — 
im[H, x\/h, where H is the Hamiltonian and m the elec- 
tron mass, one can show the equivalence of Eqs. (8a) 
and (9a) at resonance. The a; 3 dependence off reso- 



nance in length gauge can be justified by analogy with 
Abraham-Lorentz radiative damping, with a damped os- 
cillator Green function 
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Inserting the expression T n (Lo) = T n (oj) into Eq. Q, one 
obtains the length-gauge form for the imaginary part of 
the polarizability off resonance (loq — > coon, 7 — > T„). 

Quite surprisingly, the question of whether one should 
use the length or velocity forms for the decay width off 
resonance, i.e., in the interval < lo < loq h , has not 
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FIG. 1. Plot of the integrand x 5 e~ bx L s (x) defined in 
Eq. \l\h\ in the interval characteristic of the thermal peak 
x w 5/6 = 0.01583 [Fig. (a), parameter s = I] and near the 
resonant peak x ~ 0.375 [Fig. (b), indistinguishable curves 
for s — 1 and s = 3]. The thermal peak [Fig. (a)] yields 
the dominant contribution to the model integral J defined in 
Eq. dTTab. 



been answered conclusively in the literature. It has of- 
ten been stressed (e.g., in Ref. [T2]) that the electric field 
strength E is a physical observable and thus gauge in- 
variant while the gauge-dependent vector potential A is 
not. Lamb noted in footnote 88 on p. 268 of Ref. [13] that 
the interpretation of the quantum mechanical wave func- 
tion is only preserved in the length gauge with the dipole 
r ■ E interaction, because the kinetic momentum changes 
from p — > p — e A in the presence of a vector potential, 
and therefore the p operator in the quantum mechanical 
Hamiltonian of the atom cannot be interpreted any more 
as a kinetic momentum if the vector potential is nonvan- 
ishing (see also Chap. XXI of Ref. QH and Refs. QI1QS]). 
On the other hand, in Ref. [7], the authors explain in the 
text after Eq. (2) that "the velocity [gauge] form follows 
originally from the [fully relativistic] QED description" 
and should therefore be used off resonance, for the ob- 
vious reason that the nonrelativistic limit of the Dirac 
matrix vector a, which enters the relativistic expression 
for the self-energy |17) is the momentum operator p/(mc). 
While the length-gauge results seem to be generally pre- 
ferred in the literature, the use of the length versus veloc- 
ity forms remains controversial, and all numerical results 
below are therefore indicated for both velocity and length 
gauge; further considerations on the choice of the gauge 
are beyond the scope of the current article. 

Model example. — In order to illustrate the numerical 
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evaluation of Eq. ([!]) , we consider a dimensionless model 
integral which is obtained by replacing frequency, transi- 
tion width and (3 with their dimensionless equivalents: 

uj^x = hcj/E h , T^j = nv/E h , 

cj -> x = huj /E h , l3^b = f3E h = E h /(k B T), 
where Eh is the Hartree energy. The resulting integral 

J / dxx 5 e~ bx L s (x), (11a) 
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contains the "Boltzmann factor" e~ which models the 
hyperbolic sine in the denominator of the integrand of 
Eq. ([!]), and the imaginary part is taken for a single res- 
onance function that models the Lorentzian line profile 
[with s 



1 and s 



3 for the analogues of Eq. (8b) 
and Eq. (9b I respectively]. We choose the temperature 
as T — 1000 A", corresponding to b = 315.775, and res- 
onance parameters for the lowest (1S-2P) transition in 
the hydrogen atom: xq = |, and 7 = 1.5162 x 10~ 8 . The 
imaginary part of the Lorentzian "polarizability term" 
L s (x) is highly peaked near the resonance energy x — x 
and for s = 0, and the full width at half maximum 
(FHWM) is equal to 7. The relative change of the prefac- 
tor x 5 exp(-bx) over the interval x G (xq — ^ 7, xq + | 7) 
is smaller than 10~ 5 . One thus conjectures that the re- 
placement L s (x) — > ir5(x — Xq) should provide for an 
excellent numerical approximation to the contribution of 
the Lorentzian peak to the integral J. Indeed, 
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confirming that P\ « P2 for the peak term. However, 
this treatment ignores the possibility of far off-resonant 
driving of the transition for x <C xq . A numerical evalua- 
tion of the infrared thermal spectrum leads to the result 
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which is larger than ( 12 ) by roughly 30 orders of magni- 
tude. The exact numerical result for J fulfills J s sa Q s to 
six decimals for s = 1 and s — 3; the effect is dominated 
by off-resonance absorption (see also Fig. [I]). 

Simple Atoms. — Returning from our model exam- 



ple (11a) to realistic simple atoms, we provide results 
for the numerical integration of Eq. ([!]) for hydrogen 
atoms in Fig. [2] and for helium atoms in their ground and 
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FIG. 2. The characteristic slowdown time due to friction 
for ground state atomic hydrogen as a function of the black- 
body radiation temperature. The solid line shows the results 
using the dynamic polarizability of Eq. Q in length-gauge 
form [Eq . q9bf ], the long-dashed line is the velocity-gauge 
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form \8b\, the shaded area is in between, and the short-dashed 
Its from Dirac 8 peaks given in Eq. pj. 
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FIG. 3. Same as Fig. [2]for a ground state helium atom. 
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FIG. 4. Same as Figs. [2] and [|] for metastable helium (2 3 Si) 



metastable triplet states in Figs. [3] and [4j respectively. 
The friction force is expressed in terms of its correspond- 
ing characteristic slowdown time r = mv/F. For atomic 
hydrogen and helium, the dynamic polarizability Q has 
been used with the parameters listed in Tables |T] |III| 
The input data have been partially calculated by us, and 
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TABLE I. Frequencies, oscillator strengths, lifetimes, and 
level widths for transitions between the ground IS state and 
n x P levels within the hydrogen atom. 



n 


uj 0n [a.u.] 


/on [a.u.] 


l/r„ [ns] 


T n [a.u. 




2 


0.375 000 000 


0.416 196 717 


1.595 


1.5162 10" 


-8 


3 


0.444 444 444 


0.079 101562 


5.268 


4.5911 10" 


-9 


4 


0.468 750 000 


0.028 991029 


12.346 


1.950710" 


-9 


5 


0.480 000 000 


0.013 938 344 


23.949 


1.007810" 


-9 



TABLE II. Same as Table|l]for transitions between the ground 
1 S state and the singlet n P levels within the helium atom. 



If 


U0n [a.u.] 


fon [a.u.] 


l/r„ [ns] 


r n [a.u.] 


2 


0.779 881291 


0.276 1647 


0.555 


4.3514 10" B 


3 


0.848 578 015 


0.073 4349 


1.728 


1.3998 10~ 8 


4 


0.872 654 727 


0.029 8629 


3.975 


6.0852 10~ 9 
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0.883 818 387 


0.015 0393 


7.653 


3.1607 10" 9 



TABLE III. Same as Table |H| for transitions from the triplet 
metastable 2 3 S state to triplet n 3 P levels (helium atom). 



n 


W n [a.u.] 


/on [a.u.] 


l/r„ [ns] 


T„ [a.u. 




2 


0.042 065187 


0.539 0861 


98.202 


2.462310" 


iii 


3 


0.117148 294 


0.064 4612 


98.154 


2.4635 10" 


10 


4 


0.142 905 024 


0.025 7689 


142.886 


1.692310" 


10 


5 


0.154 678191 


0.012 4906 


225.643 


1.071610" 
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the transition frequencies and oscillator strengths have 
been verified against those given in Refs. [TSJ . The 
total decay rates used in the calculation include the de- 
cays to both 1 S and 1 D states. The temperature at which 
the full Lorentz profile results start to deviate from the 
Dirac-(5 peaks is given by T* = where x is the greater 
of the two real and positive (rather than complex) solu- 
tions of the equation x 7 e~ x = ffTT 5 ^ (velocity gauge) 
and x 9 e~ x — ^7r 7 ^ (length gauge, n = 2 is the princi- 
pal quantum number of the lowest excited state). For an 
equation of the form x n e~ x = A, this particular solution 
can be expressed as x = — n W-i(— A x / n /n), where W 
is the generalized Lambert W function [2U]. In velocity 
gauge, T* evaluates to 3293 K for hydrogen, 6927 K for 
singlet and 346K for triplet helium. In length gauge we 
have T* = 2954 K for hydrogen, 6208 K for singlet and 
312 K for triplet helium (confirmed in Figs. [3] and [4]). 

Conclusions. — In this Letter, we show that far off- 
resonant driving of atomic transitions yields the domi- 
nant contribution to the blackbody friction force on mov- 
ing atoms, due to the overlap of the infrared tail of the 
Lorentzian profile with the infrared thermal peak of the 
blackbody radiation. It is thus imperative to take the 
finite lifetime of the atomic resonances and their corre- 



sponding width into account. Numerical results for sim- 
ple atoms are provided in Figs. [2]— [4| The feasibility of an 
experimental verification of the predictions of this Letter 
remains to be studied. Of the atomic systems considered 
here, the largest effect is expected for the metastable 3 Si 
state in helium. In this case and for a temperature equal 
to the melting point of tungsten (3695K) the character- 
istic slowdown time is computed to be 3016s (~ 50 min- 
utes), which makes the friction effect difficult to observe 
in laboratory experiments, but perhaps not impossible. 
The general importance of an accurate understanding of 
blackbody friction for astrophysical processes has already 
been stressed in Ref . [1] . Further remarks on conceivable 
astrophysical consequences of the calculations reported 
here are beyond the scope of this Letter. 
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